Anisotropic magnetoresistance in ferromagnetic atomic-sized metal contacts 
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Recent experiments in ferromagnetic atomic-sized contacts have shown that the anisotropic mag- 
netoresistance (AMR) is greatly enhanced and has an asymmetric angular dependence as compared 
with that of bulk samples. The origin of these effects is still under debate. In this work we present 
a theoretical analysis of the AMR in atomic contacts of the 3d ferromagnetic materials. Our results 
strongly suggest that the anomalous AMR stems from the reduced symmetry of the atomic con- 
tact geometries. We also present calculations supporting the idea that the pronounced voltage- and 
temperature dependence in some experiments can be attributed to impurities near the constrictions. 
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When a metallic wire is shrunk to the atomic scale, 
its transport properties are significantly altered [1]. One 
interesting and novel example is the anisotropic magne- 
toresistance (AMR) that arises when the magnetization 
throughout a ferromagnetic device is rotated uniformly 
a certain angle 9 with respect to the current direction. 
In polycrystalline bulk samples such rotation induces a 
relative change in the conductance AG/G that varies as 
cos 2 9 with an amplitude on the order of 1% [5] . Re- 
cently Bolotin et al. [5] found that the AMR of permal- 
loy atomic-sized contacts can be considerably enhanced 
as compared with bulk samples and that it exhibits an 
angular dependence clearly deviating from the cos 2 9 law. 
Additionally, they found a significant voltage dependence 
on the scale of millivolts, which led them to interpret the 
effect as a consequence of conductance fluctuations due 
to quantum interference jj. Independently, Viret and 
coworkers [5] reported similar results in Ni contacts, but 
also the occurrence of conductance jumps upon rotation 
of the magnetization. Similar stepwise variations of the 
conductance have been found in Co nanocontacts [B]. 

The jumps have been interpreted as a manifestation of 
the so-called ballistic AMR (BAMR) [TJ. According to 
this theoretical concept the rotation of the magnetization 
in a ballistic contact could result in an additional band 
crossing the Fermi energy, leading to an abrupt change 
in the conductance on the order of e 2 /h. This has been 
theoretically shown to occur in ideal infinite chains of 
Ni and Fe [TJ. However, realistic ferromagnetic contacts 
made of transition metals are not expected to be ballistic 
[5] and so the interpretation of the conductance jumps in 
terms of BAMR is at least questionable. Indeed, Shi and 
Ralph [9] have suggested that these jumps might origi- 
nate from two- level fluctuations due to changes in atomic 
configurations jlOj . Thus, there remain important open 
questions about AMR in atomic contacts, concerning the 
origin of the enhanced amplitude, the anomalous angular 



dependence, the occurrence of conductance jumps, and 
the voltage dependence. 

In this Letter we address these questions using a combi- 
nation of a tight-binding (TB) model, molecular dynam- 
ics (MD) simulations, and a simple toy model. Our cal- 
culations suggest that the enhancement of the AMR am- 
plitude and the deviations from the cos 2 9 law in atomic 
contacts stem from spin-orbit coupling (SOC) together 
with a reduced symmetry of the contact geometry. We 
do not find signs of BAMR in realistic contact geometries. 
Finally, we find that the conductance of pure atomic con- 
tacts has no voltage dependence on the the scale of mil- 
livolts, but the addition of impurities can lead to a sig- 
nificant voltage dependence of the AMR signal. 

Our description of the AMR in ferromagnetic atomic 
contacts of the 3d metals (Ni, Co, and Fe) is based on the 
following non-orthogonal TB Hamiltonian that includes 
the SOC: 

E (^w+G^)^- (i) 

Here i,j run over the atomic sites, a, J3 denote the dif- 
ferent atomic orbitals, and cr = |, j the spin. The spin- 
polarized parameters hf} - n .„ and the spin-independent 
overlaps Siajp are taken from the parametrization of 
Refs. [TTJ [12] , which is known to accurately reproduce the 
band structure and total energy of bulk ferromagnetic 
materials and has been successfully applied to the de- 
scription of ferromagnetic atomic contacts H¥j . The 
atomic basis is formed by 9 orbitals (3d, 4s, 4p), which 
give rise to the main bands around the Fermi energy in 
Fe, Co, and Ni. The matrix elements of the intra-atomic 
SOC are given by h^°J ~, — £i{iaa\l-s\i(3a') , where land 
s are the orbital and spin angular momentum operators. 
These matrix elements depend on the spin quantization 
axis (with polar angle 9 and azimuthal angle </>), which 
we rotate uniformly to simulate the AMR experiments 
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with saturated magnetization. We use the SOC constant 
= 70 meV, which is the experimental value for Ni. 
We stress that, while this type of model is not expected 
to be as accurate as ab initio counterparts, it captures 
the essential physics. It also has the advantage that it 
allows us to analyze phenomena like BAMR in realistic 
contact geometries, which is currently outside the scope 
of ab initio methods. 

In order to describe the electronic transport we use the 
nonequilibrium Green's function formalism [15] . Briefly, 
the atomic contacts are divided into three parts, a central 
region C containing the constriction and the left/right 
(L/R) leads, which we model as infinite surfaces. The 
retarded Green's functions of the central part read 



Gcc = 



H 
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where Yl x = tcxgxx t* cx are the lead self-energies (X = 
L,R). Here, t cx = sS cx - H<&> with H^, S cx 
the hopping elements and overlaps between the C region 
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FIG. 1: (Color online) (a) Ideal Ni one-atom contact in fee 
[111] direction with atoms on lattice positions. Green atoms 
are those in the atomic constriction, yellow ones are part of 
the surfaces used to model the leads. (b,c) Channel decom- 
position and the total linear conductance as a function of 8 
for different angles 4>. The relative conductance is defined as 
G re i = G(8, 4>)/(G(8, 4>)}g — 1. (c) Channel decomposition of 
(b). (d)-(f) Same as (a)-(c), but with the contact distorted by 
randomly shifting the red atoms by up to 5% of the nearest- 
neighbor distance. 



and the lead X, and gxx is a surface Green's function. 
In general H^ ' and Sl,/?. depend on the bias voltage 
V. The ^-dependent transmission matrix is t(e, V) = 

r^ /2 G cc r^ /2 , where T x = i(E x - st^). The current 
then adopts the standard Landauer-Buttiker-likc form 



I(V) 



deT(e,V)[f L (e,V)-f R (s,V)}, (3) 



where Jl,r are the Fermi functions and r(e, V) — Tr [t^t] 
is the transmission function. The low-temperature linear 
conductance can be written as G = (e 2 /h) T m where 
T n are the transmission coefficients, i.e. the eigenvalues 
of t^t at Fermi energy ef- 

In the calculations presented here we have neglected 
the SOC in the leads. We have checked that this only 
introduces a small change in the contact resistance that 
does not alter the conclusions but reduces the computa- 
tion time enormously. On the other hand, we determine 
self-consistcntly the on-site energies of the atoms in the 
constriction by imposing the local charge neutrality that 
metallic elements should exhibit. 

We now apply this method to calculate the conduc- 
tance of Ni atomic contacts. First we consider an ideal 
geometry with the atoms kept fixed on fee lattice posi- 
tions and forming pyramid- like tips in the [111] direction 
that end in a common central atom [Fig. [it a)]. In Figs. 
[TJb,c) the conductance and its channel decomposition as 
function of 9 for several values of are shown. Surpris- 
ingly, the conductance of this one-atom contact exhibits 
the bulk-like AMR with a cos 2 8 dependence (minimum 
at 9 = 0), an amplitude of 0.5% and practically no de- 
pendence on </>. In fact, the individual channels show a 
more complicated dependence on 8 7 and the amplitude of 
variation for one channel can be bigger than that of the 
total conductance, but in the latter these features cancel 
and the cos 2 9 dependence is recovered. 

It seems obvious that the cancellation is related to the 
high symmetry of the ideal geometry. To test this idea 
we have distorted the contact by shifting randomly the 
atomic positions by up to 5% of the nearest-neighbor 
distance [Fig. [ltd)]. As seen in Fig. [TJe), the individ- 
ual channels now show roughly the same amplitude of 
variation with 9 as in the ideal contact, but due to the 
disorder they exhibit a more complex 9 dependence, and 
a strong dependence on 4>. As a consequence, the con- 
tributions of the channels no longer cancel out and the 
AMR can have a different amplitude, with the conduc- 
tance extrema shifted in 9 and with a strong dependence 
on cf> [Fig. [ljf )] . This example illustrates that the origin 
of the anomalous angular dependence and amplitude of 
the AMR in atomic contacts can be simply the reduced 
symmetry of these junctions together with the fact that 
the conductance is mainly determined by a few atoms 
in the narrowest part of the constrictions. We want to 
point out that we have found similar results for Co and 
Fe atomic contacts, which confirms this conclusion |15j . 
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FIG. 2: (Color online) Contact evolution of a Ni junction grown in the fee [001] direction as obtained from the simulations 
of Ref. [14]. (a) Spin-projected, Gf,|, and total conductance in the absence of SOC and total conductance averaged over 
6,cj> in the presence of SOC. Vertical lines correspond to the contact geometries in (b). Inset: relative AMR amplitude 
AG/{G)e = (G max ,e(0) — Gmin,e(4>))/{G(8, 4>)}e vs. inverse averaged conductance, (c) Conductance vs. 6 for the geometries in 
(b) and with <f> in steps of tt/6. (d) Relative AMR [G(6, <f>)/{G(0, <t>))e,4> - 1] in % on a "Bloch sphere", (e) Same as (b)-(d) for 
a thick regular geometry with 324 atoms with the relative AMR values multiplied by 10 for visibility. 



Since the geometry plays such a prominent role in the 
AMR, it is important to determine the geometries that 
can be realized in an actual experiment. For this pur- 
pose, we have carried out classical MD simulations of the 
formation of Ni atomic contacts, following Ref. [14]. An 
example of the contact evolution is shown in Fig. [2] Here, 
we start with an ideal Ni bar containing 112 atoms on 
lattice sites in fee [001] direction. The bar is attached to 
rigid surfaces that are separated in a stepwise manner, 
simulating the elongation process of a break junction. 
In Fig. |2ja) we show the evolution of the spin-projected 
and total conductance during elongation in the absence 
of SOC. Adding SOC introduces only a small change in 
the averaged total conductance. As usual, the sudden 
atomic rearrangements are reflected as steps in the con- 
ductance p] . The vertical lines refer to the geometries of 
Fig. [2jb) obtained during the elongation. For them we 
have computed the dependence of the conductance on 9 
and (f>, and the results are shown in Fig. [2^c). For con- 
tact 1, which is just an elastic deformation of the ideal 
contact, the conductance has two types of behavior de- 
pending on 4>: one is cos 2 0-like, while the other is clearly 
more complex. In order to visualize the overall angu- 
lar dependence, we show the relative AMR on a "Bloch 
sphere" [Fig. [2^d)]. The contact 1 and hence its AMR 
have an approximate four-fold symmetry. When defor- 
mations emerge in the contact, the angular dependence 
becomes irregular and strongly dependent on (j). For ex- 
ample, for contact 2 there is a strong variation of AMR 
with </>, and depending on its value, the AMR amplitude 
can be almost one order of magnitude larger than in the 
bulk limit or cancel almost entirely. As the contact evolu- 



tion proceeds, the AMR has an amplitude of around 2%. 
In almost all cases, the conductance is not only shifted in 
9, but it also has a more complicated behavior than just 
cos 2 9. Approaching the tunnel regime (G < 0.1e 2 //i), 
we do not observe a further increase of the AMR ampli- 
tude, contrary to experiments [3]. One reason may be 
that the isolated tip atoms in tunneling regime exhibit 
a finite orbital moment [16 (not considered here), which 
may lead to a local deviation of the spin-quantization 
axis from the field direction and an additional increase 
in resistance. Finally, in the limit of thick contacts, we 
recover bulk behavior with an amplitude of 0.45%, as 
shown in Fig. [2je). The statistical analysis of all contact 
geometries shows an increase of AMR to 2% on average in 
the last steps before breaking, see inset of Fig|2^a). This 
confirms that the lack of symmetry in atomic contacts 
gives rise to the enhancement of the AMR signal. 

We have not found indications of tip resonances, which 
are present in ideal one-dimensional geometries |17l 118] 
and which were suggested as the origin of the experimen- 
tal findings [17]. For example, the transmission for the 
contact of Fig. [TJa) has almost no structure around ep 
on the scale of millielectronvolts, as shown in Fig. |3ja). 
We thus believe that the voltage and temperature depen- 
dences reported in Refs. [5J [TO] are indeed associated with 
impurities close to the constriction, as reported earlier for 
non-magnetic junctions [HE [20]. Following Ref. [2U], one 
may estimate that the reported [3j voltage period of a 
few millivolts can stem from impurities located hundreds 
of nanometers away from the contact. Such length scales 
cannot be modeled realistically, but we have developed a 
toy model to support the idea. As represented schemat- 
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FIG. 3: (Color online) (a) r(e, V = 0) for the contact of 
Fig. [TJa) with 6 1 = 0,tt/2. (b) Sketch of the chain model, (c) 
corresponding transmission for 6 — 0,7r/2 in the absence of 
impurities, (d) r(e, V = 0) with an impurity N = 751 sites 
from the scattering region, and (e) relative nonlinear conduc- 
tance Grei(6,V) = G{6,V)/{G(e,V)}e-l- (f) Relative linear 
conductance vs. 6 and (g) voltage dependence of nonlinear 
conductance at 6 — it/2 for indicated temperatures. Results 
at K in (f,g) correspond to the dashed lines in (e). 



ically in Fig. [3jb) , we model the contact with two semi- 
infinite linear chains connected at the tips, and describe 
the system with a nearest-neighbor TB Hamiltonian. For 
simplicity we only include the d xy and d x i_ y i orbitals 
(the E 2 doublet), with a spin-splitting of A = 0.5 eV. In 
the absence of SOC there are four independent equiva- 
lent chains with onsite energies ±A/2, hopping integrals 
tti = 1.0 eV in the leads, and tt c = 0.4 eV between the 
tip atoms. The SOC that couples the chains is restricted 
to the two tip atoms, with a coupling constant £ = 0.2 
eV. For further details, see Ref. [T5]. As we show in Fig. 
[3jc), there is no structure in r(e, V = 0) and the AMR 
is small due the high symmetry of the geometry. In this 
example ep — 0.42 eV, which makes the sign of the AMR 
opposite to that in Fig. [3|a). We note that the results of 
Fig. [TJb) are qualitatively reproduced [15] , 

Things change drastically when an impurity is intro- 
duced in one of the chains. We model it by using 
slightly reduced hoppings n s — 0.8 eV between two atoms 
N = 751 sites away from the constriction in the right 
lead. Figure (3^d) shows the resulting transmission func- 
tions for 6 = 0,7r/2. In addition to the expected oscil- 
lations in energy, there appears a modulation depending 
on 9 that stems from the interference with the impurity. 
When a finite voltage V is applied over the tip atoms, 
the nonlinear conductance G(V) = dl/dV exhibits os- 
cillations as a function of both voltage and the angle 8. 
The relative variations shown in Fig. [3^e) bear a striking 



resemblance to the experimental results of Ref. |3J. Fi- 
nally, we show in Figs.[3jf,g) the temperature dependence 
of both the linear and non- linear conductance. The effect 
of temperature is to smooth the K characteristics, again 
compatible with the experimental observations |10j . 

In summary we have shown that the anomalous mag- 
nitude and angular dependence of the AMR in ferromag- 
netic atomic-sized contacts can be explained naturally in 
terms of the reduced symmetry of the atomic junction 
geometries. We predict a strong anisotropy of the con- 
duction channels, but have not found any signature of 
BAMR. We have presented a simple model which illus- 
trates that the pronounced voltage and temperature de- 
pendence found in some experiments may originate from 
the presence of impurities close to the constriction. 
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